Excitons in carbon nanotubes: an ab initio symmetry-based approach 
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The optical absorption spectrum of the carbon (4, 2) nanotube is computed using an ab initio 
many-body approach which takes into account excitonic effects. We develop a new method involving 
a local basis set which is symmetric with respect to the screw symmetry of the tube. Such a 
method has the advantages of scaling faster than plane-wave methods and allowing for a precise 
determination of the symmetry character of the single particle states, two-particle excitations, and 
selection rules. The binding energy of the lowest, optically active states is approximately 0.8 eV. 
The corresponding exciton wavefunctions are delocalized along the circumference of the tube and 
localized in the direction of the tube axis. 
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A wealth of extraordinary results concerning the me- 
chanical and electrical properties of carbon nanotubes 
(NTs) have been reported in the last few years Q . Until 
recently, however, their optical properties have not re- 
ceived the same attention: experimental work has often 
been hindered by low emission efficiency, and the inter- 
pretation has been complicated by the fact that tubes of 
different species and orientation are normally mixed to- 
gether in the same sample, making it difficult to assign 
the measured spectra to a single NT species. 

Very recent experiments have indicated that these lim- 
itations can be overcome. Improved optical efficiency 
has been obtained by isolating NTs in porous materi- 
als 0], in solution Q, or on patterned substrates Q. It 
was possible to assign optical spectra to specific NTs via 
their characteristic vibrations in resonant Raman [{| or in 
near-field experiments |fj ■ The observation of electrically 
induced optical emission from a carbon NT FET Q has 
paved the way for a new class of single-molecule exper- 
iments and devices. These advances therefore establish 
optical spectroscopies as powerful characterization tools 
for NTs. Nanotubes also hold great promise for novel 
nanoscale opto-electronic and photonic applications 
because the optical gap of NTs spans a very large fre- 
quency range which overlaps with the range of interest 
in the field of telecommunications. 

In spite of such fervent interest in this subject, the fun- 
damental nature of optical excitations of NTs is not yet 
understood. The possible relevance of excitonic effects in 
these systems was pointed out in a pioneering paper by 
Ando [8j. In general, it is well known that the electron- 
hole interaction plays a crucial role in one-dimensional 
systems, not only in the ideal case 0, but also in realis- 
tic systems such as semiconductor quantum wires jl(l or 
polymer chains 1 1 1 1] where excitons dominate the optical 
spectra. The binding energies are however very sensitive 
to the spatial extent of the single-particle wavefunctions 
and to the (anisotropic) dielectric screening [J^. In the 
case of NTs, one might expect these quantities to be sen- 



sitive to size and geometry. Moreover, given the peculiar 
nature of the electronic states, even the smaller diameter 
NTs cannot be regarded as pure one-dimensional sys- 
tems. It is therefore extremely difficult to estimate the 
actual relevance of excitonic effects without a realistic 
calculation of their optical excitations. 

In this Letter, we present an efficient ab initio approach 
to compute the spectra, exciton energies and wavefunc- 
tions of carbon NTs. Such an approach includes the full 
three-dimensional dependence of the electron states and 
interactions. Our formalism exploits a fully symmetrized 
Gaussian basis set and allows us not only to reduce the 
system size and computation time considerably, but also 
to profit from the selection rules involving a new quantum 
number which comes into play in the symmetry charac- 
terization of the problem. 

As a preliminary step, a calculation based on the 
density-functional theory within the local density ap- 
proximation jl.'{l| is performed with a plane-wave basis, 
pseudopotentials and supercells and the band states 
"0nfc(G) are calculated in Fourier space. These states are 
then projected onto a set of basis functions which are 
symmetrized sums of Gaussian orbitals centered on the 
atoms. The basis is constructed to be simultaneous eigen- 
states of the commuting operators G\ and G2 ■ G± repre- 
sents a discrete translation of the tube in the z-direction 
by the length of the cell T. G2 is the screw-symmetry 
operation consisting of a combined rotation and trans- 
lation: a rotation by an angle of 2ir/N, where N is the 
number of hexagons in the supercell; and a translation 
along the z-direction by (M/N)T, where M depends on 
the geometry and size of the tube (for details see Ref. Q]). 
The basis functions have the form 
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The functions {</> n /m(r)} are Gaussians defined as 
n *m(r) = r l exp(-r 2 /al)Yi m (r). They have been used 
with considerable success in similar calculations for ex- 
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tended systems |l5|. N is the number of hexagons per 
unit cell, a — 1,2 labels one of the two inequiva- 
lent atoms in the unit cell. / q = e 27riq and = 
e 2-ni(Mq+h)/N are fa e eigenvalues of the operators G\ and 
G2, respectively, h is a new quantum number, analogous 
to the azimuthal quantum number m, which ranges from 
1 to N. 

The quantities needed for the calculation of both 
the random phase approximation (RPA) dielectric func- 
tion lla and the Bethe Salpeter equation (BSE) ma- 
trix Il7i fisl are the Gaussian transition elements 



j^nk,n'k' = (nk\x~ k ' :hd \n'k'). They are given by 
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where ip" k is the Gaussian component of the wave- 
functions of the system [Greek letters a, (3, 7, . . . denote 
the quartuple (a,n,l,m)], and = h(nk) — h(n'k') is 
the difference between the quantum number h for the 
states nk and n'k 1 . The three-basis-function integrals 
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are given by 



FIG. 1: Band structure in the extended zone scheme. Each 
panel represents a full Brillouin zone (FBZ). The solid line 
represents a continuation of a band with character ft = 1 in 
the first Brillouin zone, and the dashed line begins with a 
band with character h — 2. The top of the valence band is 
represent by a horizontal dotted-dashed line. 
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The exchange term is: 



The phase-multiplied three-center integrals £™V 1 
calculated from the pure three-center integrals I^^g 1 
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1 " are the basis vectors (atomic positions in the 
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reduced cell) for Gaussians a, (3, and 7. 

The evaluation of the Gaussian transition elements and 
the three-basis-function integrals in Eq. (J2J) and Eq. © 
is the most expensive part of the calculation and, for a 
given precision, scales like ~ (NNk) 2 ~ R 2 , where Nk is 
the number of fc-points and R is the radius of the tube. 
This scaling is considerably faster than that of a plane- 
wave calculation, which scales like ~ TR 4 \og(TR 2 ). Our 
method has the additional advantage of scaling indepen- 
dently of the cell length T and allows for a calculation of 
tubes like the (4,2) NT, which have unusually long cell 
lengths. 

The expressions for the BSE matrix are simple func- 
tions of the matrix elements in Eq. J5J). We distinguish 
two cases: (a) If the light is polarized along the tube axis, 
z, the selection rule for h is h = h! . For this polarization 
we include only those transitions which conserve h. In 
so doing, the BSE matrix is reduced to N independent 
blocks, (b) The same holds for circularly-polarized light 
with polarization vector e± = l/v / 2(e a: ± ie y ) with the 
only difference that the selection rule is h = h! ± 1. The 
direct term for both polarizations is: 
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where h(ck) — h(vk) = (±1) for the e z (e±) polariza- 



tion, and are the screened and bare Coulomb 

interactions expressed in the symmetrized basis of Gaus- 
sians (the Einstein sum convention is used to sum over 
the Greek indices) calculated following the procedure in 
Ref. The screening is calculated in RPA. Building 

the matrices in Eq. Q and Eq. (J5J does not require much 
computation time because the Greek indices range from 
1 to the number of Gaussians (fifty in this work) which 
is small and does not depend on system size. 
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The exciton wavefunction for the 



excitation is: 
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where the expansion coefficients ^"1 satisfy the BSE: 
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^vk,c'v'k' ~ * y cvk,c'v'k'- Since we 
are interested only in the effect of the electron-hole 
interaction on the spectra, we use LDA energies for 
the E n k of Eq. without quasiparticle corrections [2(| . 
The imaginary part of the dielectric function is given 
by e 2 (w) = 4ir/iu 2 J2n\ M n\ 2 $(u - n„) where M n = 

J2cvk^cvk( v M v ' e |ck). In this work, the momentum 
operator is used in lieu of the velocity operator since the 
former is easier to express in Gaussians. The difference 
in another carbon-based system, i.e., diamond, has been 
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FIG. 2: The left panels show the RJDOS and the optical ab- 
sorption for the e 2 polarization and the right panels for the e+ 
polarization. In both bottom panels, the solid line represents 
the spectrum with electron-hole interaction and the light grey 
line the spectrum in the single-particle picture. The dotted- 
dashed line represents the spectrum with exchange term only. 
All spectra are computed with a broadening of 0.07 eV, with- 
out self-energy corrections. 

shown to be a 15 % to 20 % effect in the peak intensity 
of absorption spectrum [2l]], but has a negligible effect 
on the peak position. 

We focus on the carbon (4,2) NT 22]. Fig. □ shows 
its electronic band structure. It is represented in the 
extended zone scheme to clearly show how the bands un- 
fold and how they are related to the quantum numbers 
h and k. Every panel corresponds to a FBZ. The entire 
band structure is shown with light grey lines. The thick 
solid (dashed) lines are lines of constant character h = 1 
(h = 2). In this scheme, it is clear how to represent the 
set of all optically permitted transitions for both polar- 
izations (e z and e±). For the e 2 polarization, only the 
vertical transitions from the solid valence line to the solid 
conduction line and those from the dashed valence line 
to the dashed conduction line are optically allowed. For 
the e± polarization, only the vertical transitions from the 
solid to dashed line, or viceversa, depending on the he- 
licity of the polarization, are optically allowed. Due to 
symmetry, both helicities ( e + and e_) yield the same 
spectrum. In the non-interacting theory, these transi- 
tions can be represented using a reduced joint density 
of states (RJDOS). It is obtained by excluding from the 
joint density of states all optical transitions forbidden by 
symmetry. The RJDOS for both polarizations is shown 
in the top two panels of Fig. [5J The corresponding ab- 
sorption spectra are shown in the bottom two panels. 
They are computed with 24 k-points in the FBZ and a 
broadening of 0.07 eV. The BSE matrix computed with 
this k-point sampling and for the energy range of inter- 



FIG. 3: The top panel shows a portion of the band structure 
in the extended zone scheme with the relevant bands [solid 
bold line (h — 1) in Fig. 0. The bottom three panels show 
the weighted contributions of the transitions which comprise 
three excitonic states that correspond to the first three peaks 
in the absorption spectrum for the e z polarization. 



est is small with dimension w 300x300. For the e z po- 
larization, the absorption spectrum with excitonic effects 
(solid line of bottom left panel of Fig.0) has three peaks: 
the lowest two peaks correspond to bound excitons with 
binding energies Eb = 0.8 eV and Eb = 0.2 eV, while the 
third one is an unbound resonance. This is consistent 
with the value of 1.0 eV found for an (8,0) nanotube in 
another work [23|. 

To gain a better understanding of the nature of the 
excitonic states giving rise to these three peaks, we con- 
sider the three corresponding excitonic states with the 
highest oscillator strength M n . Fig. shows, in the ex- 
tended zone scheme, the weighted contributions of the 
optical transitions to these states (bottom three panels). 
State 1 is made up of transition A shown in the band di- 
agram (top panel) . State 2 is derived from transitions A 
and C, and state 3 from transitions A and B. Transitions 
A, B, and C correspond to the Van Hove singularities in 
the RJDOS labeled in the top left panel of Fig. H Note 
that the resonance represented by state 3 is coupled with 
transitions in the continuum region, appearing in Fig. [3| 
(bottom panel) as two sharp peaks each corresponding 
to a single Bloch vector. 

From the bottom left panel of Fig. one can directly 
compare the absorption with both direct and exchange 
terms included in the BSE (solid line), the absorption 
with exchange term only (dotted-dashed lines), and the 
absorption without excitonic effects (the shaded line): it 
is evident that excitonic effects radically alter the absorp- 
tion spectrum. Moreover, most of the effect comes from 
the direct term, as the exchange term alters the spectrum 
by a small amount. 
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FIG. 4: Exciton wavefunctions for the two lowest, optically 
active excitons of the (4,2) carbon NT. The top panel shows 
the Ah = exciton (Eb = 0.8 eV) and the bottom panel 
the Ah — 1 exciton (Eb = 0.6 eV). The color plots represent 
the projection on the tube lateral surface of the probability 
of finding the hole when the electron is fixed in the origin 
slightly above one of the atoms, z is the tube axis. 9 is the 
circumference direction. 



In the bottom right panel of Fig. [3 we represent the ab- 
sorption spectrum for the e + polarization: the spectrum 
calculated using the non-interacting theory is greatly 
suppressed when both exchange and direct terms are in- 
cluded in the BSE Hamiltonian. This phenomenon was 
discovered in another theoretical work on nanotubes \2m 
using time-dependent local density approximation |25| . 
To see if such a suppression is to be attributed to the 
exchange term [equivalent to local field effects (LFE)], or 
to the direct electron-hole Coulomb attraction, we plot 
an additional curve with the exchange term only (dotted- 
dashed line). Fig. [3 shows that the suppression is largely 
due to the exchange term and that the direct term affects 
the spectrum less, in agreement with Ref. |24J. The main 
effect of the direct term is to open up a significant peak 
below the lowest van Hove singularity, corresponding to a 
bound, Ah = I exciton with binding energy E b — 0.6 eV. 

We plot the exciton wavefunction 0W(r e ,rj,) [see 
Eq. |JB}] by fixing the position of the electron: Fig. 0] 
shows the two lowest lying, optically active excitons for 
the e z (Ah — 0) and the e + polarization (Ah — I). Both 
excitons are localized in the direction of the tube axis, z, 
but delocalized along the tube circumference, 9. This 
behavior is in agreement with the variational calculation 
of Ref. j2(| . The wavefunction for the even parity (and 
hence, optically inactive), Ah = exciton located 0.1 eV 
below the first absorption peak (not shown) is similar to 
the lowest optically active (odd parity) exciton. 

In conclusion, we have developed an efficient method 
for performing BSE calculations on NTs of all sizes and 
chiralities that uses a local, symmetry-based approach. 
Applying this method to a (4,2) NT, we have shown 
that excitonic effects are important and have quantita- 
tively determined these effects. For the e 2 polarization 
(Ah = 0) , the excitonic effects give rise to three enhanced 



peaks, and the lowest excitation is 0.8 eV below the onset 
of the single-particle continuum. The corresponding ex- 
citon is localized along e z , but delocalized along the tube 
circumference. For the e + polarization(A/i = 1), the ex- 
change part of the BSE matrix drastically suppresses the 
absorption spectrum. The direct part opens up a large 
peak 0.6 eV below the single-particle singularity corre- 
sponding to a bound, localized exciton. 
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